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ABSTRACT 


An  algorithm  Is  presented  which  computes  optimal  weights  for  arbitrary 
linear  arrays.  The  application  of  this  algorithm  to  In  sUu  optimal  reshading 
of  arrays  with  failed  elements  Is  discussed.  It  Is  shown”  that  optimal 
reshading  can  often  regain  the  original  sjdelobe  level  by  slightly  Increasing 
the  malnlobe  beamwldth.  Three  examples  are  presented  to  Illustrate  the 
algorithm's  eTfectlveness.  Hardware  and  software  Issues  are  discussed. 
Execution  time  for  a  25  element  array  Is  typically  between  one  and  two  minutes 
on  a  HP9836C  microcomputer. 
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I.  INTRODUCTION 

A  linear  array  of  discrete  elements  (sensors)  often  experiences  element 
failures  in  situ.  These  failures  can  significantly  Increase  the  sidelobe 
levels  of  the  array  wavenumber  response,  depending  on  how  many  elements  fail 
and  where  the  elements  are  located  within  the  array.  We  discuss  here  an 
optimal  reshading  (reweighting)  algorithm  which  can  be  applied  in  situ  to 
reduce  the  sidelobe  levels  to  the  original  design  level.  In  many  common 
element  failure  situations,  optimal  reshading  can  regain  the  original  sidelobe 
level  by  slightly  increasing  the  malnlobe  beamwidth.  In  arrays  which 
experience  significant  element  failures,  optimal  reshading  is  still  possible, 
but  may  be  of  limited  use.  Three  examples  given  below  demonstrate  a  few  of 
the  possibilities. 

An  algorithm  for  optimal  reshading  was  first  proposed  in  [1]  by  Strelt 
and  Nuttall.  Their  algorithm  utilized  the  general  purpose  subroutine  [2]  to 
solve  a  specially  structured  "linear  programming"  problem.  Unfortunately, 
their  algorithm  required  hours  of  computation  time  and  large  amounts  of 
computer  storage  on  a  minicomputer  (the  VAX  11/780)  to  optimally  reshade  a  50 
element  array  with  five  failed  elements.  Consequently,  their  algorithm  is  not 
useful  for  in  situ  optimal  reshading. 

The  shading  algorithm  proposed  here  differs  from  Strelt  and  Nuttall's 
primarily  in  that  we  solve  their  linear  programming  problem  using  a  new 
general  purpose  subroutine  [3,4],  herein  referred  to  as  Algorithm  635. 
Algorithm  635  uses  the  special  structure  of  the  linear  programming  problem  to 
reduce  time  and  storage  requirements  by  orders  of  magnitude.  Algorithm  635 
can  be  Incorporated  easily  In  Strelt  and  Nuttall's  original  approach.  A 
significant  algorithmic  improvement  was  discovered  in  the  course  of  this  study 
and  is  described  below.  The  resulting  shading  algorithm  is  fast  enough  and 
small  enough  to  execute  successfully  on  microcomputers  (such  as  the  HP9836C 
used  here)  in  only  a  few  minutes.  Typical  execution  time  for  a  25  element 
array  is  under  two  minutes;  for  a  50  element  array,  execution  time  is 
typically  under  10  minutes.  The  current  algorithm,  and  the  HP9836C  with  its 
Inherent  transportability,  comprise  an  effective  system  for  optimal  reshading 
in  situ. 


II.  OPTIMAL  ARRAY  SHADING 

The  wavenumber  response  of  a  linear  array  composed  of  N  discrete 
omnidirectional  elements  located  at  arbitrary  fixed  positions,  xn,  is  given 
by 


T(k)  =  2^  wn  exp[-ikxn]  (1) 

n=l 


where  wn  are  the  element  weights  and  the  independent  variable  k  denotes 
wavenumber  in  radians  per  unit  length.  The  element  weights  are  required  to  be 
real,  but  this  entails  no  loss  of  generality  (see  below  in  section  III). 

Also,  from  (1),  T(-k)=T*(k)  for  real  weights  (asterisk  denotes  conjugation). 


1 


TM  No.  861004 


so  it  is  unnecessary  to  consider  negative  values  for  k  and  we  confine  our 
attention  to  non-negative  k. 

The  array  response  as  a  function  of  k  can  be  considered  to  be  composed  of 
a  mainlobe  beamwidth  and  a  sidelobe  region.  The  objective  of  the  optimization 
process  is  to  make  |T(k)|  as  small  as  possible  on  the  user  specified  sidelobe 
Interval.  Array  weights  which  achieve  this  objective  are  said  to  be  optimal. 
The  optimization  process  usually  produces  equi-valued  sidelobes  in  the 
sidelobe  region. 

Weights  that  are  optimal  for  a  full  array  do  not  remain  optimal  after  the 
array  experiences  element  failures.  To  partially  compensate  for  failed 
elements,  the  array  Is  optimally  reshaded  by  undertaking  the  optimization 
process  again  and  incorporating  knowledge  of  which  elements  have  failed.  As 
the  examples  below  will  show,  the  effectiveness  of  this  strategy  depends  upon 
how  many  elements  have  failed  and  the  location  of  these  elements  in  the  array. 

The  sidelobe  Interval  Is  defined  differently  depending  on  the 
Inter-element  spacing  of  the  array.  For  an  array  with  periodically  spaced 
elements  and  no  failures,  the  sidelobe  interval  is  defined  to  be 
[Kq,(2w/D)-Kq],  where  Kg  Is  calculated  from  the  desired  sidelobe  level 
and  the  number  N  of  array  elements.*  0  Is  the  physical  distance  from  sensor 
to  sensor.  Furthermore,  the  minimization  Interval  can  be  reduced  to 
[Kq,*/D],  since  the  response  of  this  array  Is  symmetric  about  k=*/D.  Kq 
Is  typically  the  point  on  the  mainlobe  response  which  Is  equal  in  magnitude  to 
that  of  the  sidelobes,  but  this  Is  not  always  true  for  seriously  degraded 
and/or  aperiodic  arrays  (see  example  3  below).  For  arrays  with  aperlodically 
spaced  elements,  the  sidelobe  Interval,  denoted  by  [Kg.K^],  must  be  chosen 
by  Inspection  of  a  nonoptimal  beampattern  or  some  other  means.  |T(k)|  must  be 
minimized  over  the  full  [Kg.K-i]  range  since,  in  general,  an  aperiodic 
array  response  Is  not  symmetric  about  any  wavenumber  other  than  k=0.  The 
ability  to  specify  arbitrary  Kg  and  K-j  Is  particularly  useful  for  those 
applications  involving  aperlodically  spaced  elements  because  lower  sidelobe 
levels  may  be  obtained  by  looking  at  different  minimization  regions. 

The  optimization  process  deals  with  element  failures  In  an  array  In  the 
following  way: 

Step  1.  Maintain  mainlobe  beamwidth  and  permit  the  sidelobe  levels  to 
rise. 

Step  2.  Regain,  If  possible,  the  original  sidelobe  level  by  broadening 
the  mainlobe. 

Broadening  the  mainlobe  by  increasing  Kg  (step  2)  Is  performed  only  If  the 
sidelobe  level,  even  after  optimal  reshading,  has  risen  to  an  unacceptable 
value  because  of  element  failures.  Thus,  step  1  is  normal  algorithmic 


*  For  an  N-element  array  and  -t  dB  peak  sidelobes,  we  have  K(y=(2/D)arccos(l/Zn) 
where  2Zg=[r+(r2-l)1/2]'/'M  +  [r-(r2-l  )1/2]1/M,  r=10t/2°,  and  M=N-1 . 

The  Inter-element  spacing  0  Is  assumed  to  be  half  of  the  so-called  design 
wavenumber,  and  N  Is  the  number  of  array  elements  before  failures. 
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procedure,  and  step  2  requires  some  iteration  In  specifying  Kq  and/or  K-j 
because  a  compromise  has  to  be  made  between  the  mainlobe  beamwidth  and  the 
level  of  the  sidelobes. 

The  solution  of  the  array  problem  In  the  original  formulation  [1]  is 
mathematically  equivalent  to  solving  an  overdetermined  system  of  complex 
linear  equations.  Unacceptably  high  sidelobes  result  If  this  system  is  solved 
In  the  usual  least  squares  sense,  so  It  Is  necessary  to  solve  the  system  so 
that  the  magnitude  of  the  maximum  residual  error  is  minimized.  There  now 
exists  [3]  an  efficient  algorithm  and  corresponding  FORTRAN  code  [4]  for 
solving  problems  of  this  sort  to  high  accuracy. 

To  obtain  the  beamformer  equation  in  an  appropriate  format  to  utilize 
this  algorithm,  we  normalize  the  peak  response  of  T(k)  so  that  T(0)=1.  This 
gives 


£  “n  '  1  '  <2> 

n=l 


We  solve  <q.  (2)  for  the  Nth  weight,  w^,  and  substitute  in  eq.  (1)  to  obtain: 


N-l 

T(k)  =  exp[-ikxN]  +  £  wn[exp(-1kxn)  -  exp(-ikxN)]  (3) 

n=l 


By  sampling  T(k)  at  the  M  equi-spaced  points 


[K  -  V 

km  =  Kq  +  —  (m-1)  ,  m=l , . . . ,  M , 


(4) 


we  can  write  the  problem  of  minimizing  the  peak  sidelobe  level  of  the  array 
response  as 


min  max 
jwl  l<m<M 


N-l 


f™-I 


a«m  Wn 

mn  n 


n=l 


(5) 


where  the  complex  numbers  fm  and  amn  are  defined  by 


f  =  exp[-1kmx.,3 
m  m  NJ 


amn  =  exp[-1kmxN]  '  exp["ikmxn] 


(6) 
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The  problem  (5)  is  precisely  the  form  necessary  for  application  of  Algorithm 
635.  For  theoretical  details  of  this  algorithm,  the  Interested  reader  is 
referred  to  [3]. 

Sometimes  a  few  of  the  optimum  weights  for  arrays  with  failed  elements 
are  observed  to  be  negative,  particularly  those  on  the  end  elements.  If  the 
weights  are  applied  in  hardware,  providing  a  180  degree  phase  factor  on  the 
element  output  may  not  be  desirable  or  possible.  However,  Algorithm  635 
allows  the  selection  of  all  non-negative  weights;  this  Is  Implemented  by  the 
addition  of  constraints  to  eq.  (5).  Usually,  but  not  always,  an  element  Is 
zeroed  If  It  would  have  had  a  negative  weight.  From  eq.  (2)  It  follows  that, 
If  all  the  element  weight  values  are  required  to  be  positive,  they  must  be 

between  0  and  1.  The  requirement  that  weights  w-j . wN_^  be  between  0 

and  1  can  be  written  mathematically  as 


,  n=l . N-l 


(7) 


Algorithm  635  requires  these  N-l  constraints.  Algorithm  635  can  also 
Incorporate  any  number  of  general  constraints  of  the  form 


S  bmn  wn  "  cm 
n 


where  cm  and  dm  are  constants.  The  requirement  that  w^  also  be 
non-negative  gives 


<  d 
-  m 


(8) 


(1 


N-l 

£  wn>  -  7 
n=l 


or 


N-l 

£ 


n=l 


w_ 


(9) 


which  Is  clearly  a  special  case  of  the  general  constraints. 


III.  ALGORITHM  IMPROVEMENTS 

Several  changes  to  the  algorithm  presented  In  [1]  enable  significant 
reduction  In  the  need  for  computational  intensity.  Lewis  and  Streit  [5]  have 
proved  that,  for  a  general  line  array  shaded  so  that  It  has  optimal  sldelobe 
levels  when  steered  through  the  same  number  of  degrees  either  side  of 
broadside,  there  exists  a  set  of  optimal  weights  which  are  real.  Thus, 
complex  weights  do  not  need  to  be  considered.  This  fact  allows  an  approximate 
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eight-fold  reduction  in  computation  time  and  a  two-fold  reduction  in  storage 
requl rements. 

It  Is  clear  that  the  50  element  example  run  in  Streit  and  Nuttall  [1]  was 
significantly  oversampled  In  wavenumber.  Their  beampattern  can  be  reproduced 
with  a  four-fold  reduction  In  the  sampling  of  T(k)  (see  Example  2  below),  and 
this  in  no  way  detracts  from  the  practical  application  of  the  algorithm.  A 
significant  reduction  in  computation  time  is  realized  by  decreasing  the 
number,  M,  of  beampattern  samples. 

A  significant  algorithmic  modification  made  to  Algorithm  635  further 
decreases  computation  time.  We  have  labeled  this  modification  "fast  costing" 
and  it  is  an  Important  step  in  making  the  algorithm  feasible  on  microcomputers 
such  as  the  HP9836C.  In  order  to  properly  describe  this  modification,  some 
familiarity  with  the  simplex  method  of  linear  programming  and  reference  [3]  is 
assumed. 

Algorithm  635  can  be  broken  into  two  fundamental  computational  operations 
called  "costing"  and  "pivoting."  "Costing"  determines  the  so-called  "minimum 
reduced  cost  coefficient"  and  requires  2NM  multiplies,  where  N  is  the  number 
of  discrete  array  elements  and  M  Is  the  number  of  samples  taken  of  the 
beampattern.  "Pivoting"  is  a  basis  update  and  requires  N2  real  multiplies. 

It  is  clear  that  the  speed  of  the  algorithm  is  intimately  related  to  the 
number,  M,  of  samples  taken  of  the  beampattern,  as  well  as  the  number,  N,  of 
discrete  array  elements.  Since  M  is  larger  than  N,  "costing"  requires  more 
multiplies  than  "pivoting." 

"Costing"  In  the  linear  array  application  means  that,  in  each  simplex 
Iteration,  the  "discretized  absolute  value"  of  every  sidelobe  sample  of  the 
wavenumber  response  function,  TJkm),  m=l ,  ...,  M,  is  computed  to  determine 
the  "minimum  reduced  cost  coefficient"  of  the  current  "basic  feasible 
solution."  By  proceeding  through  a  finite  sequence  of  such  "basic  feasible 
solutions,"  we  arrive  at  the  solution  of  the  "discretized  problem."  As  shown 
In  [3],  this  implies  that  the  computed  optimal  wavenumber  response  function 
can  have  sidelobe  levels  that  are  theoretically  at  most  .04  dB  higher  than  the 
true  optimum  sidelobe  level.*  "Fast  costing"  refers  simply  to  the  fact  that 

we  first  determine  which  of  the  sidelobe  samples  T ( kq,) ,  m=l . M,  has  the 

largest  true  absolute  value,  and  then  compute  the  "discretized  absolute  value" 
of  this  one  complex  number.  Therefore,  only  one  "discretized  absolute  value" 
calculation  is  performed  in  each  simplex  iteration  instead  of  M  such 
calculations.  The  resulting  reduction  in  computational  effort  Is  significant 
in  microcomputing  environments.  The  drawback  is  that  the  use  of  "fast 
costing"  prevents  the  simplex  algorithm  from  converging  to  a  solution  of  the 
"discretized  problem."  Fortunately,  however,  it  can  be  proved  that  we  must 
approximate  the  solution  in  a  well  defined  sense.  In  the  linear  array 
application,  "fast  costing"  results  in  the  computed  optimum  beampattern  having 
sidelobe  levels  that  are  theoretically  at  most  .08  dB  higher  than  the  true 


*  The  theoretical  error  of  at  most  .04  dB  is  derived  by  taking 
20  log-|o(sec(ir/p)) ,  where  p=32.  The  term  sec(ir/p)  is  the  error  bound 
discussed  in  [3]. 
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optimum  level.*  This  is  a  small  price  to  pay  for  major  execution  time 
improvements. 


IV.  ALGORITHM  IMPLEMENTATION  FOR  IN  SITU  USE 

An  algorithm  must  be  reliable,  easy  to  use,  and  fast  when  executing  on 
portable  microcomputers,  to  be  useful  for  In  situ  application.  The  following 
section  details  the  most  important  hardware  and  software  Issues  addressed  to 
enable  in  situ  optimal  reshading  of  arrays  with  failed  elements. 

The  algorithm  has  been  coded  In  BASIC  and  is  comprised  of  Algorithm  635 
and  an  array  processing  driver  program.  Algorithm  635  solves  the  linear 
program  for  a  set  of  optimal  weights,  given  data  supplied  by  the  driver 
program.  The  driver  performs  the  initial  setup  based  on  several  user  Inputs 
and  provides  all  program  output. 

The  driver  program  may  be  used  with  linear  arrays  having  either  periodic 
or  aperlodically  spaced  elements.  Program  output  consists  of  a  graph  of  the 
optimal  beampattern,  a  graph  of  the  optimal  normalized  element  weights,  and 
several  parameters  pertinent  to  the  specific  problem.  Provision  is  made  for 
storing  the  weights  in  a  separate  data  file  for  possible  use  with  digital 
beamformers. 

An  HP  specific  software  modification  was  made  by  setting  up  the  Input 
data  arrays  (equation  6)  In  buffers  so  that  they  are  accessible  for  a 
one-dimensional  multiply.  For  large  array  dimensions,  Indexing  a  doubly 
subscripted  data  array  and  performing  a  dot  product  takes  more  time  on  the 
HP9836C  than  reading  In  a  data  array  from  a  buffer,  doing  a  MAT  multiply,  and 
performing  a  summation.  (A  MAT  multiply  Is  simply  an  element  by  element 
multiply  of  two  equally  dimensioned  data  arrays).  However,  this  procedure  Is 
more  time  consuming  when  the  Input  data  arrays  are  very  small  (l.e.,  the 
number  of  elements  In  the  line  array  Is  small).  The  breakeven  point  occurs  at 
around  12  or  13  elements,  so  it  was  decided  to  incorporate  this  speed 
enhancement  for  the  longer  running,  larger  line  arrays  and  tradeoff  some  speed 
reduction  on  the  smaller  line  arrays. 

To  obtain  fast  execution  times  for  in  situ  application,  we  use  one 
hardware  speed  enhancement,  a  12.5  MHz  fast  CPU  card  with  16  kbytes  of  cache 
memory.  This  hardware  supplement  is  available  from  HP  for  use  on  the 
HP9836C .  Cache  memory  is  fast  memory  resident  on  the  CPU  card  for  quick 
Instruction  acquisition.  The  use  of  the  fast  CPU  board  rather  than  the  8  MHz 
clock  present  In  the  standard  computer  configuration  results  In  an  approximate 
factor  of  two  Increase  in  observed  speed. 

The  complete  program  is  precompiled  by  use  of  software  and  a  floating 
point  math  card  available  from  the  INFGTEK  company.  Precompllation  reduces 
most  computational  portions  of  the  BASIC  code  to  machine  language,  giving  an 
additional  three-fold  reduction  In  computation  time.  It  is  also  desirable  to 
upgrade  the  operating  system  for  the  HP  to  its  latest  revision.  All  work  on 


*  Fast  costing  squares  the  error  bound,  giving  sec?(ir/p),  or  .08  dB  when 
p=32. 
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these  problems  was  run  using  the  BASIC  3.0  operating  system  and  the  hardware 
supplements  noted  above. 

Computation  time  is  defined  as  time  spent  in  Algorithm  635  and  does  not 
Include  the  small  amount  of  setup  time  required  by  the  driver  program. 
Computation  times  are  for  the  compiled  BASIC  program  run  on  the  HP9836C  with 
the  special  hardware  additions  mentioned  above. 

The  program  described  here  needs  just  over  303  kbytes  of  Internal  memory 
in  addition  to  the  memory  required  by  the  operating  system  to  execute  on  the 
HP9836C .  This  is  the  amount  of  space  required  by  fixing  the  maximum  array 
size  at  N=50,  and  allowing  at  most  M=256  beampattern  samples.  Users  can 
change  dimensions  to  suit  their  specific  needs,  but  storage  requirements 
presently  are  directly  proportional  to  the  product  NM.  Even  for  a  much  larger 
number  of  line  array  elements,  it  is  unlikely  that  memory  restrictions  would 
prove  to  be  a  problem  on  the  HP9836C  since  extra  memory  boards  of  1  Mbyte  each 
are  readily  available. 

Ongoing  modifications  should  further  enhance  the  capability  and  speed  of 
the  BASIC  algorithm  and  driver.  The  addition  of  the  ability  to  handle 
directional  sensors  is  both  useful  and  straightforward  to  implement. 

Execution  of  Identical  code  on  the  new  HP  300  series  computers,  which  have  a 
16.6  MHz  clock  rate,  should  further  reduce  the  computation  time.  Computation 
times  on  the  order  of  5  minutes  for  a  50  element  array  and  one  minute  for  a  25 
element  array  are  anticipated. 

It  Is  possible  to  run  the  BASIC  program  In  Its  uncompiled  state.  The 
execution  of  the  program  with  cache  memory  and  the  fast  CPU  board  as  the  only 
enhancements  results  in  computation  times  of  approximately  25  minutes  for  a  50 
element  array  and  4.5  minutes  for  a  25  element  array. 

A  copy  of  the  entire  program  Is  available  from  the  authors.  Our  specific 
implementation  in  HP  BASIC  utilizes  several  hardware  and  software  devices  to 
achieve  computational  efficiency,  some  of  which  may  not  be  pertinent  to  other 
BASIC  operating  systems  running  on  comparable  machines.  Users  will 
undoubtedly  find  It  necessary  to  make  modifications  to  the  code  to  allow  it  to 
run  on  other  HP  equipment  or  In  BASIC  on  the  VAX. 


V.  EXAMPLES 

The  following  examples  demonstrate  the  utility  of  the  current  algorithm 
for  application  in  situ  and  provide  Insight  into  different  situations  that 
might  arise  when  reshading  equispaced  arrays  with  failed  elements.  If  optimal 
reshading  can  restore  the  array's  original  design  sidelobe  level  by  slightly 
increasing  the  mainlobe  beamwldth,  than  we  say  that  the  optimal  reshading  has 
been  effective.  Optimal  reshading  Is  effective  in  many  common  element  failure 
situations.  When  the  array  is  severely  degraded,  optimal  reshading  is  less 
effective  but  is  still  useful  in  reducing  the  negative  Impact  of  element 
failures.  These  examples  demonstrate  that  the  effectiveness  of  reshading 


7 


TM  No.  861004 


depends  upon  the  number  of  element  failures  as  well  as  the  location  of  the 
failed  elements  within  the  array. 

Missing  elements  are  modeled  by  zeroing  the  appropriate  weights.  In 
these  examples,  N  refers  to  the  number  of  intact  array  elements,  M  Is  the 
number  of  beampattern  samples,  and  Kq  Is  calculated  by  using  the  equation  in 
an  earlier  footnote.  We  define  the  malnlobe  width  to  be  twice  Kq  In  all 
three  examples. 


EXAMPLE  1.  Effective  reshading 

This  example  demonstrates  that  reshading  can  restore  the  original 
sldelobe  level  of  an  array  response  by  slightly  Increasing  the  malnlobe 
beamwldth.  In  a  25  element  equispaced  array,  originally  designed  for  -30  dB 
sldelobes,  elements  2  and  4  have  failed.  Therefore,  N=23,  M=128,  and 
Kq=.6877.  We  first  keep  the  malnlobe  width  fixed  and  allow  the  sldelobe 
level  to  rise.  See  Figure  1.  The  peak  sldelobe  level  has  risen  to  -26.86  dB 
below  the  malnlobe  and  the  malnlobe  width  Is  unchanged.  If  the  sldelobe  level 
after  reshading  Is  too  high,  an  alternative  to  discarding  or  repairing  the 
array  Is  to  broaden  the  malnlobe  beamwidth.  In  Figure  2,  Kq  Is  Increased  to 
.775  and  the  peak  sldelobe  level  diminishes  to  -30.04  dB  below  the  malnlobe. 

A  tradeoff  must  always  be  made  between  an  enlarged  malnlobe  beamwidth  and  an 
acceptable  peak  sldelobe  level.  In  this  case  the  malnlobe  was  Increased  12.7 
percent  In  order  to  recover  the  original  sldelobe  level.  Execution  times  on 
the  HP9836C  are  between  one  and  two  minutes  for  Figures  1  and  2. 


EXAMPLE  2.  Moderately  effective  reshading 

This  example  Is  taken  from  Strelt  and  Nuttall  [1].  Because  of  the 
Improvements  detailed  In  Section  III  above,  the  current  algorithm  runs  faster 
on  the  HP9836C  than  on  the  VAX  11/780,  although  the  floating  point  multiply 
time  on  the  HP  In  Its  basic  configuration  Is  roughly  200  times  slower  than  on 
the  VAX. 

Consider  a  linear  array  with  50  equispaced  elements,  Initially  designed 
for  peak  sldelobes  of  -30  dB  relative  to  the  malnlobe.  Figure  3  shows  the 
classical  Dolph-Chebyshev  beampattern  with  -30  dB  sldelobes  throughout  the 
minimization  range  C Kq . ( 2ir/0 ) -Kg ] .  This  was  computed  using  the  current 
algorithm  In  6.11  minutes.  (This  ideal  case  could  have  been  computed 
analytically.) 

Now  we  suppose  that  five  elements,  7,22,40,43,50,  of  the  array  have 
failed.  The  optimal  response  after  reshading  the  array  Is  shown  In  Figure  4. 
The  peak  sldelobe  level  has  risen  to  -25.51  dB,  but  we  have  maintained 
malnlobe  beamwidth  and  retained  full  steering  capability.  In  this  example 
N=45  and  M=128. 

This  example  (Figure  4)  took  7.47  minutes  on  the  HP9836C  and  required  292 
simplex  Iterations.  The  algorithm  of  Strelt  and  Nuttall  required  38.4  minutes 
and  402  iterations  on  the  VAX. 
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The  recovery  of  the  original  sidelobe  level  is  possible  (Figure  5).  The 
mainlobe  beamwidth  must  be  Increased  by  the  large  factor  257.6  percent 
(Kq=.871)  and  the  execution  of  this  task  takes  8.98  minutes  and  requires  351 
Iterations.  The  constraint  that  all  the  weights  lie  between  0  and  1  Is  used. 
It  Is  necessary  to  use  the  constraint  In  this  instance  because  otherwise  a 
dislocation  of  the  maximum  response  from  k=0  results.  This  dislocation  is  due 
to  the  presence  of  too  many  negatively  weighted  elements. 


EXAMPLE  3.  A  severely  degraded  array 

This  example  shows  that,  for  severely  degraded  arrays,  recovery  of  the 
original  sidelobe  level  may  not  be  possible  by  Increasing  the  mainlobe 
beamwidth,  even  after  optimal  reshading.  Consequently,  control  of  the  level 
of  the  first  sidelobe  must  be  relinquished  In  order  to  gain  control  of  the 
level  of  the  remaining  sidelobes. 

Consider  a  25  element  array  with  elements  11  and  14  failed.  The  original 
sidelobe  level  Is  -30  dB.  Here  N=23,  M=128,  and  Kg=.6877.  Figure  6  shows 
the  algorithm's  optimal  response  to  this  configuration.  It  is  a  significant 
observation  that,  in  this  case,  small  perturbations  of  Kg  will  not  affect 
the  level  of  the  sidelobes.  Only  when  the  first  sidelobe  is  Incorporated  Into 
the  mainlobe  beamwidth  (Kg=1.27)  does  the  level  of  the  remaining  sidelobes 
return  to  the  original  desired  value  (see  Figure  7).  It  Is  apparent  that 
decreasing  the  minimization  interval  by  moving  Kg  far  enough  to  the  right 
will  Improve  the  approximation,  but  one  must  give  up  control  of  the  first 
sidelobe  to  reduce  the  others  to  acceptable  levels.  The  net  effect  of  losing 
two  elements  so  close  to  the  center  is  that  negligible  emphasis  Is  placed  on 
the  remaining  center  elements  (12  and  13)  and  the  rest  of  the  aperture  Is 
reshaded  as  If  it  were  two  separate  arrays. 

This  situation  cannot  be  overcome  by  using  different  weights.  The 
optimal  property  of  the  array  problem  formulation  and  solution  tells  us  that 
no  weights  exist  which  can  suppress  all  the  sidelobes  below  a  certain  level. 
Thus,  this  array  has  lost  too  many  elements  and  performance  cannot  be  restored 
to  its  original  design  levels  merely  by  reshading. 

We  have  chosen  to  relinquish  control  of  the  first  sidelobe  to  gain 
control  of  the  level  of  the  remaining  sidelobes.  We  pick  the  first  sidelobe 
merely  for  ease  of  Implementation;  modification  of  the  algorithm  to  forfeit 
control  of  a  different  sidelobe  could  also  have  been  done.  The  need  to 
relinquish  control  of  the  first  sidelobe  level  has  only  appeared  in  cases  of 
severe  array  degradation  due  to  element  losses. 


VI.  CONCLUSIONS 

Arrays  that  have  failed  elements  can  be  reshaded  to  obtain  optimal  array 
response  functions.  Optimal  reshading  Is  effective  In  many  common  element 
failure  situations.  When  the  array  Is  severely  degraded,  reshading  Is  less 
effective,  but  still  can  be  used  to  reduce  the  negative  Impact  of  element 

failures. 
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Optimal  reshading  can  be  accomplished  In  situ,  quickly  and  reliably,  on 
portable  microcomputers  using  the  algorithm  described  here.  Arrays  with  25 
elements  routinely  run  In  less  than  two  minutes  and  computation  time  for  a  50 
element  array  Is  less  than  10  minutes.  The  algorithm  can  be  applied  to  arrays 
of  evenly  or  unevenly  spaced  linear  geometry. 

The  above  examples  (and  others)  support  the  generally  accepted  notion 
that  failure  of  near-center  elements  Is  more  detrimental  to  the  array  response 
than  failure  of  near-edge  elements. 


Another  application  of 
three  dimensional  geometry, 
probably  will  depend  upon  N 
samples)  In  the  same  manner 


Algorithm  635  Is  to  arrays  of  planar  and  arbitrary 
Computation  times  for  these  more  general  arrays 
(number  of  sensors)  and  M  (number  of  beampattern 
as  for  linear  arrays. 
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Figure  1.  Optimized  array  response  and  normalized  weights 
for  25  elements  with  elements  2  and  4  missing. 
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Figure  2. 


Array  response  and  normalized  weights  for  Example  1  with  Ko=.775. 


13 


TM  No.  861004 


K0= .3381  K  (RRDXMETER)  Kl=12.2283 


Figure  3.  Classical  Dolph-Chebyshev  array  response  and  normalized  weights 

for  N--50  and  -30  dB  sidelobes. 
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Figure  4.  Optimized  array  response  and  normalized  weights  for 
50  elements  with  elements  7,  22,  40,  43,  and  50  failed. 
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Figure  5.  Recovery  of  original  sidelobe  level,  Example  2,  Kq=.871. 
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Figure  6.  Optimal  array  response  and  normalized  weights 
for  25  elements  with  elements  11  and  14  failed. 
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Figure  7.  Recovery  of  original  sidelobe  level,  Example  3  with  Kq-1.27. 
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